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We study the temperature dependence of the strange and charm quark chemical 
potentials in the phase diagram of nuclear matter, within a modified and gener¬ 
alized hadron gas model, in order to consider phase transitions and to describe 
phenomena taking place outside the hadronic phase. We employ, in a phenomeno¬ 
logical way, the Polyakov loop and scalar quark condensate order parameters, 
mass/temperature-scaled partition functions and enforce flavor conservation. We 
propose that the resulting variation of the heavy quark chemical potentials can be 
directly related to the quark deconfinement and chiral phase transitions. Then, 
the chemical potential of the strange and charm quark can be considered as an 
experimentally accessible ’’order parameter”, probing the phase diagram of QCD. 


1 Introduction 

One of the main problems in the study of the phase transitions, occurring on the 
level of strong interactions, is finding an unambiguous observable, which would act 
as an experimentally accessible ’’order parameter” [1]. All proposed QGP signa¬ 
tures (strangeness enhancement, J/V’ suppression, dileptons, resonance shift and 
broadening, etc.) have already been observed in heavy ion collisions, however, we 
have seen, that they also occur, to some extent, m p — p oi p — A interactions 
where no QGP production is theoretically expected. The physical quantity needed 
should exhibit a uniform behavior within each phase, but should change when a 
critical point is reached in a phase transition. It has been earlier suggested [2-4] 
that the chemical potential of strange quarks may be the sought-for macroscopic 
and therefore measurable thermodynamic quantity. The case of [2-1-1] flavors was 
thoroughly studied and it was shown that the change in the sign of the strange 
quark chemical potential, from positive in the hadronic phase to negative in the 
deconflned phase, may indeed be a unique indication of the deconfinement phase 
transition. Here we will review the basic aspects of the model and present the [2-1-2] 
flavors version, which is a generalization of the model with the inclusion of c-quark 
and charm hadrons. 


2 Hadronic phase 

Assuming that the system has attained thermal and chemical equilibration of four 
quark flavors (u, d, s, c), the partition function for the hadronic gas is written in the 
Boltzmann approximation: 

InZHciT, y A„ Ae) = InZ^G + (1) 
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where 


ZnZ“’^ = Z„ + Z„(A3 + A-3) (2) 

= ZK{XqK^ + A-'A,) + Zy{\1\s + A-"A;1) 

+^H(AgAs + A^ ^Ag ^) + Zn(As + A^ (3) 

and 

= Zd{XcX-^ + A,A-i) + ZdAXcX^^ + A’^A,) 

+^A<,,i:e(AcAq + A^ ^Ag + Z'B.AXqXsXc + A^ ^A^ ^A^ 

+-^ne(AsAc + Ag ^A^ (4) 

is the partition function for the non strange, strange and charm sectors, respectively. 
The charm sector also includes strange/charm mesons and baryons that lead to a 
coupling of the fugacities Ac, A^. For simplicity we have assumed isospin symmetry 
Xu = Xd = Xq, while the one particle Boltzmann partition function is given by: 

3 

The summation in Eq.(||) runs over the resonances of each hadron species with 
mass TOj, and the degeneracy factor gj counts the spin and isospin degrees of 
freedom of the j-resonance. For the strange hadron sector, kaons with masses up 
to 2045 MeV/c^, hyperons up to 2350 MeV/c^ and cascades up to 2025 MeV/c^ 
are included, as well as the 12“ states at 1672 MeV/c^ and 2252 MeV/c^. For the 
charm hadron sector, we include purely charm mesons D'^,D~,D^ and baryons 
(Ac, Ec) as well as strange-charm mesons (Df) and baryons (Sc, He) which contain 
both heavy quark flavors. All known charm resonances are taken into account with 
masses up to 2.7 GeV/c^. To derive the Equation of State (EOS) of the hadron gas 
phase we simultaneously impose flavor conservation, 

<N,-Nj> = ^^InZHoiT, V, Xq, A„ Ac) = 0 (6) 

V UfAs 

< A^c - A% > = L^lnZHciT, V, Xq, A«, Ac) = 0 (7) 

V UfAc 

which reduce to a set of coupled equations: 

Zk{X-^Xs - XqX-A + ZriXlXs - Xq^X^A + 2^h(A,A2 - A'^A^^) 
+3ZniXl - XjA + ZsAXqXAc - X-^Xj^X^A + 2ZnAX% - Xj^X^A = 0 

( 8 ) 

ZoiXaXq^ - XqX-A + ZdAXcX-^ - X-^Xs) + Za^,yAXcXI - A-'A-^) 

+^h.(A,A«Ac - A-iA.-^A-i) + Za,(A^c - A.-^a-I) = 0 (9) 

The above conditions, define the relation between all quark fugacities and tem¬ 
perature in the equilibrated primordial state. In the HG phase with finite net 
baryon number density, the chemical potentials fiq, g,s and gc are coupled through 


cf2002proc: submitted to World Scientific on February 1, 2008 2 




the production of strange and charm hadrons. Due to this coupling > 0 in 

the hadronic domain. A more elegant formalism describing the HG phase is the 
Strangeness-including Statistical Bootstrap model (SSBM) [5,6]. It includes the 
hadronic interactions through the mass spectrum of all hadron species, in contrast 
to other ideal hadron gas formalisms. The SSBM is applicable only within the 
hadronic phase, defining the limits of this phase. In the 3-flavor case, the hadronic 
boundary is given by the projection on the 2-dimensional (T, fig) phase diagram 
of the intersection of the 3-dimensional bootstrap surface with the strangeness- 
neutrality surface (/is = 0). Note that the vanishing of /ig on the HG borderline 
does not apriori suggest that /is = 0 everywhere beyond. It only states that the 
condition /is = 0 characterizes the end of the hadronic phase. Figure ^ exhibits the 
hadronic boundary for two heavy quark flavors, obtained by imposing the condi¬ 
tions /is = 0 and /ic = 0 to Eq’s. (^), (9). Observe, that there exists an intersection 
point, at Tint ^ 130 MeV and /i“* ^ 325 MeV. For an equilibrated primordial 
state (EPS) above this temperature, i.e T > 130 MeV, and low pg values, we ob¬ 
serve that as the temperature decreases, the condition /if, = 0 is realized before the 
vanishing of /tg (case I), whereas for T < 130 MeV and high pg, the opposite effect 
takes place (case II). This behavior, may be of some importance towards a possible 
experimental identification of a color superconducting phase, which is realized at a 
low temperature and high density region of the phase diagram (case II). 



Light quark chemical potential (MeV) 


Figure 1. The critical curves /ig = 0 and = 0. We distinguish two cases depending on the 
location of the equilibrated primordial state (EPS). 
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3 Chirally symmetric QGP phase 


The partition function for a four flavor Quark Gluon plasma has the form, 
. „ \ G [37 . 9„9 . Mo 


lnZQGp{T, V, y + 


+ + ( 10 ) 


where m®, is the current strange and charm quark masses respectively. Flavor 
conservation within the QGP phase yields Ag = Ac = 1 or 

^QGP(T, /Xg, /Xc) = fig, fis) = 0 (11) 

throughout this phase. Here the two order parameters, the Polyakov loop < L > 
and the scalar quark density < x/u/i > , have reached their asymptotic values. Note 
that the chirally symmetric quark-gluon plasma phase always corresponds to a 
vanishing heavy quark chemical potential. 


4 Deconfined Quark Matter phase of [2+2] flavors 

We argue that, beyond the hadronic phase, an intermediate domain of deconfined 
yet massive and correlated quarks arises, according to the following qualitative pic¬ 
ture: The thermally and chemically equilibrated primordial state at finite baryon 
number density, consists of the deconfined valance quarks of the participant nucle¬ 
ons, as well as of <7 — g pairs, created by quark and gluon interactions. Beyond but 
near the HG boundary, T > T^, the correlation-interaction between q — q is near 
maximum, ««(+) <1, a prelude to confinement into hadrons upon hadronization. 
With increasing temperature, the correlation of the deconfined quarks gradually 
weakens, ois{T) 0, as color mobility increases. The mass of all (anti)quarks de¬ 
pends on the temperature and scales according to a prescribed way. The initially 
constituent mass decreases with increasing T > Td, and as the DQM region goes 
asymptotically into the chirally symmetric QGP phase, as T —s- T^, quarks attain 
current mass. Thus, we expect the equation of state in the intermediate DQM 
region to lead to the EoS of the hadronic phase, Eq. (1), at T < T^, and to the 
EoS of the QGP, Eq. (6), at T ~ In order to construct an empirical parti¬ 
tion function for the desciption of the DQM phase, we use (a) the Polyakov loop 

< L >~ = Rd{T > Td) = 0 —> 1 as T=Td and (b) the scalar density 

< >= RxiT > Td) = 1 ^ 0 as T=Td T^. The first describes the quark 
deconfinement while the latter is associated with the quark mass scaling.We assume 
that above the deconfinement temperature, quarks retain some degree of correla¬ 
tion and can be considered as hadron-like states. Therefore, near Td a hadronic 
formalism may still be applicable. This correlation/interaction gradually weakens, 
as a result of the progressive increase of color mobility. Each quark mass scales, 
decreasing from the constituent value to the current one as we reach the chiral 
symmetry restoration temperature (T ^ T^). Thus, we consider a temperature 
dependent mass for each quark flavor, approximated by: 

m}{T) = Rx(T){mf - m°f) + m°f (12) 
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where to/ and to° are the constituent and current quark masses respectively (the 
values TO° = bMeV.mP^ = 9MeV,m° = nOMeV,m° = l.lGey have been used). 
In the same spirit, we approximate the effective hadron-like mass: 

rriiiT) = RxiT){mi-m°)+m° (13) 

where rrii is the mass of each hadron in the hadronic phase and to? is equal to 
the sum of the hadron’s quarks current mass (for example m!^ = 175MeV,m~ = 
350MeV). In the partition function of the DQM phase, the former scaling is em¬ 
ployed through the mass-scaled QGP partition function InZgQp, where all quark 
mass terms are given by Eq. , while the latter is used in the mass-scaled hadronic 
partition function InZ^Q, where all hadron mass terms are given by Eq.([l^). Em¬ 
ploying the described dynamics, we construct an empirical partition function for 
the DQM phase. 


lnZDQM{V,T,{\j}) = [I- Rd{T)]lnZ*Ha{V,T,{\s}) 

+ Rd{T)lnZ*QGp{V,T,{\f}) (/ = g,s,c) 

The factor [1 — Rd{T)] describes the weakening of the interaction of the decon- 
fined quarks, while the factor Rd{T) can be associated with the increase of color 
mobility as we approach the chirally symmetric QGP phase. The DQM partition 
function is a linear combination of the HG and QGP mass-scaled partition func¬ 
tions together with the general demand to describe both confinement and chiral 
symmetry restoration asymptotically. Note that below the deconfinement critical 
point T < Td, Rd{T) = 0, leading to InZ^QM = InZnc (with constituent quarks), 
whereas at the chiral symmetry restoration temperature T ^ T^, Rd(T) = 1 and 
InZjyQM = InZQGP (with current quark masses). In order to acquire the EoS of 
the DQM phase, we impose again the strangeness and charm neutrality conditions, 
leading to the set of equations respectively. 


and 


[1 - RdiT)] [z^(A.A-i - a,a;1) + z;.iXAl - a;'a-2) 
+2Zi{XlX, - X-^X-^) + 3Z*^{Xl - X-^) + Zi^{X,XsX, - X;^X-^X-^) 

+2Z(^^(A2Ae-A;2A-i)] +R4T)gsmfK2[^){Xs-X-^) = 0 (14) 


[1 - Rd{T)] [ZU^A-^ - + ^d.(AcA;^ - A-^AQ + 

- a - ia -2) + ZZSX,XsX, - x-^x-^x-^) 
+Z*^S>^lX, - A7"A-1)] + Rd{T)g,mfK^{^Yx, - A'^) = 0 (15) 


which must be solved simultaneously. Note that because of the strange/charm 
hadrons Ds, Sc, Dc there exists a coupling between the heavy quark fugacities As, Ac. 
By solving the above equations, for a given chemical potential g,q, we derive the 
variation of the strange and charm quark chemical potentials with temperature in 
the phase diagram. 
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5 Results for finite chemical potential 


In the case of 3-flavors and finite density, we had neglected all terms involving c- 
quarks (Ac = !)■ In this case, only the variation of fj,s was considered and Figure || 
was derived. We observe that the strange quark chemical potential attains positive 
values in the hadronic phase, becomes zero upon deconfinement, it grows strongly 
negative in the DQM domain and finally returns to zero as the QGP phase is ap¬ 
proached. It is important that fXg behaves differently in each phase, as this is what 
we are looking for from the beginning in the search for an experimentally accessible 
” order parameter”. The change in the sign of fig from positive in the hadronic phase 
to negative in the deconfined is an unambiguous indication of the quark deconfine¬ 
ment phase transition, as it is independent of assumptions regarding interaction 
mechanisms. In the case of [2-1-2] flavors the situation is slightly modified. Figure 
H exhibits the variation of the two correlated heavy quark chemical potentials with 
the temperature of the primordial state, as given by Eq’s (|^, (|l^). We observe 
that both are initially positive and then grow negative, although the change in their 
sign is realized at different temperatures, for example = 0 at TJ ~ 190MeV, 
while /Tc = 0 at T| ~ 215 MeV for a fixed value of the fugacity Xq =0.48. However, 
this difference can be easily understood if we consider Figure ^ in the framework of 
Figure As already discussed in Sec. ||, for an equilibrated primordial state (EPS) 
with T > Tint and sufficiently low /ig, /ic becomes zero earlier than as the sys¬ 
tem approaches hadronization (see Figure This is the reason why > TJ in 



Figure 2. Variation of with the temperature in the case of [2+1] quark flavors and different 
approximations or parameterizations of the order parameter R^{T). 


Figure ||. For sufficiently high /r, values and low temperatures, the opposite effect 
is present, i.e /Xc changes it’s sign at a lower temperature than the strange quark 
chemical potential. The magnitude of the difference [TJ — TJ’j, will depend on the 
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exact location of the state in the phase diagram. The fact that the fig, vanish at 
different temperatures, at the end of the respective hadronic domain, has further 
consequences, as it implies that there exists a quark ’’deconfinement region” rather 
than a certain critical line. 



Figure 3. Plot of the strange and charm quark chemical potentials in the phase diagram for 
Xq = 0.48. Notice that the change in their sign is realized at different temperatures. 


6 Experimental data 

Over the last years, data from several nucleus-nucleus collisions have been analyzed 
within thermal statistical models, employing the canonical and grand-canonical 
formalisms [7-11]. Table 1 summarizes some of the results for the quantities T, 
/ig and fig, which have been deduced after performing a fit to the experimental 
data. Figure ^ shows the phase diagram with the Ideal Hadron Gas (IHG) and 
SSBM fj,s = 0 lines, as well as the location of the mean (T, /r,) values obtained for 
every collision. We observe that all interactions studied, are consistently situated 
inside the hadronic phase, defined by the IHG model and exhibit positive fig. The 
sulfur-induced interactions, however are situated slightly beyond the hadronic phase 
defined by the SSBM. IHG calculations exhibit deviations from the SSBM as we 
approach the critical deconfinement point T = ~ 175 MeV, where the S-S and 

S-Ag interactions are roughly located. Within the IHG model the condition fig = 0 
is satisfied at a higher temperature (T ~ 200 MeV), extending the hadronic phase 
to a larger region as can be seen in Figure |[ As a consequence, fig changes sign at 
a higher temperature also and this is the reason why p,s > 0 in the analysis of [11], 
although a temperature above deconfinement (according to the SSBM) is obtained. 
Therefore, an adjustment of the IHG curve to the SSBM boundary and a new fit 
to the data are necessary [12]. The data from RHIG at -^5=130, 200 AGeV are not 
included in our discussion, since at such high energies fiq is very small and fig ^ 0 
throughout the phase diagram. The observation of negative heavy quark chemical 
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potential requires a finite baryon density system. 



Light quark chemical potential (MeV) 


Figure 4. (T, /iq) values of several interactions and their location in the phase diagram. The lines 
correspond to the hadronic boundary within the SSB and IHG models. 


7 Conclusions 

On the basis of the present analysis, we conclude that the heavy quark chemical 
potentials behave differently in each region (HG-DQM-QGP) of the phase diagram 
and, therefore, they can serve as a probe of the phase transitions. This is the first 
proposal of such an experimentally accessible ’’order parameter” that holds for a 
finite baryon density state. The appearance of negative values of /is and /tc, is a well- 
defined indication of the quark deconfinement phase transition, at T=Td, which is 
free of ambiguities related to microscopic effects of the interactions. It is important 
to add, that the observation of negative heavy quark chemical potentials would 
be also a clear evidence for the existence of the proposed DQM phase, meaning 
that chiral symmetry and deconfinement are apart at finite density. Until now, 
there is no known argument from QGD that the two transitions actually occur 
at the same temperature. Au-l-Au collisions at intermediate energies, for example 
30 < y/s < 90 AGeV, should be performed to experimentally test our proposals. 
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